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TECHNICAL NOTE D-76O 

SKID LANDINGS OF AIRPLANES ON ROCKER-TYPE FUSELAGES 
By Wilbur L. Mayo 


SUMMARY 


A study is made of the landing of an airplane on a fuselage with 
"planned" curvature of its lower surface. Initial contact is con- 
sidered to stop the vertical motion of a point remote from the center 
of gravity, thus causing rocking on the curved lower surface which con- 
verts sinking-speed energy into angular energy in pitch for dissipation 
by damping forces. Analysis is made of loads and motions for a given 
fuselage shape, and the contours required to give desired load histories 
are determined. Most of the calculations involve initial contact at the 
tail, but there are two cases of unflared landings with initial contact 
at the nose. The calculations are checked experimentally for the tail- 
low case. 

Even for design sinking speeds the calculated accelerations are 
frequently less than 1 g; thus skid landings on specially shaped fuse- 
lages appear to be quite feasible and merit consideration for special 
applications. In the calculations for tail-low landings the sinking 
speed ranges from 3 to 29 ft/sec, the maximum vertical acceleration from 
0.25g to 10. Ug, the maximum nose-down angular acceleration from 0.25 to 
17.I radians /sec^, and the maximum nose-up angular acceleration from 
0.43 to 21 radians/sec^ . In the nose-low cases a sinking speed of 
9 ft/sec gives maximum ground forces 15 and 20 percent of the airplane 
weight. Experimental check of the calculations and observations on 
lateral behavior not predicted by the calculations were obtained in 
tail- low landings of a l/lO-scale dynamic model. In cross winds repre- 
senting up to 60 mph full scale, the maximum lateral deviation of the 
model represented about 60 feet full scale. 


INTRODUCTION 


The first airplane to achieve powered flight made skid landings on 
its bottom. In the development of aviation there have been recurrent 
instances of airplanes, after employing special means for take-off, 
landing on skids. With the development of rocket take-off, vertical 
risers, and boost-glide vehicles there appeared to be ample reason for 
reinvestigation of the merits of landing on skids . Such investigation 
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appeared to be particularly desirable if the lover surface of the air- 
plane could be used so that the skid would not be required to be a special 
device mounted on a shock-absorbing strut. It has appeared that this 
might be accomplished, without requiring a deformable landing surface, by 
suitable sloping of the lower surface of the fuselage so as to convert 
sinking- speed energy into angular energy in pitch for dissipation largely 
by aerodynamic damping. The investigation of this possibility forms the 
basis of the present report. The unexplained statement in reference 1 
that the long curved skid is undesirable, and the fact that skid design 
appears to have been an art rather than a science, added to the challenge 
of the investigation. Aside from possible application of its results to 
special-purpose airplanes, it was thought that a contribution might be 
made to the problems of the wheels-up landing and the hard landing in 
which collapse of the landing gears is accompanied by impact of the 
fuselage . 

The investigation was hampered by difficulty in assuming undersurface 
shapes which would give desired load variations, and difficulty in assuming 
landing conditions and load histories for which reasonable shapes of the 
fuselage undersurface would be calculated. These difficulties were largely 
overcome by fixing variables at contact of the center of the fuselage and 
then employing both forward and backward integration for desired load vari- 
ations. This procedure had the disadvantage that both the undersurface 
shape and the landing conditions were results of the calculation, and 
thus ordinary methods of analysis in which comparisons are made for the 
same initial conditions or, with other conditions held constant, a single 
variable is changed, could not be used. Therefore an increased number of 
calculations and a more complicated analysis of the effects of major vari- 
ables were necessary. Readers with limited interest in the subject may 
omit this analysis and may otherwise select data and draw their own con- 
clusions. For this purpose the data summary given in table I, arranged in 
the same order as the discussion and figures, may be useful. These data 
are applicable to airplanes which are tail- low at landing speeds. Other 
data, more limited in scope, may be applied to airplanes which normally 
land with their fuselage level, if they execute nose- low unflared con- 
tacts in order to rock on the curved fuselage undersurface. 


SYMBOLS 


^l'CpjC^C^ coefficients in differential equations, defined by 
equations (26) to (29) 
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aerodynamic pitching-moment coefficient about airplane 
center of gravity, CMa/pY^Sc 
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mean aerodynamic chord, ft 

air drag on airplane, lb 

energy due to vertical motion of weight not supported by 
wing lift 

force; ground force normal to the plane of the ground 
when used without subscript, lb 

vertical force applied by external device such as landing 
gear or arresting gear, lb 

acceleration due to gravity, ft/ sec^ 

air lift on airplane, lb 

horizontal distance from contact point to center of 
gravity, ft 

length from rear of fuselage undersurface to center of 
gravity, measured along reference line of the fuselage, 
ft 

moment, lb- ft 

moment applied by external device such as landing gear 
or arresting gear, lb-ft 




mass of airplane , slugs 

ratio of Fg to the airplane inertia, Fg^my 
radius of curvature, ft 

radius of curvature of the rocker at the center position 
beneath the center of gravity, ft 

radius of evolute circle from which a thread might be 
unwound to construct the involute rocker on each side 
of the center of gravity, ft 

radius of curvature of the rocker at the tail point, ft 

radius of gyration in pitch of airplane, ft 
wing area, sq ft 
generali zed trim, r\i 


time, sec 

landing speed, ft/sec 
weight of airplane, lb 

distance from center of gravity to point of contact of 
rocker with ground, measured along reference line of 
fuselage, positive for rearward contact points, ft 
(fig. l) 



generalized distance of center of gravity from ground, 

yu/r 


height of center of gravity above ground, ft 
(fig. 1) 


distance from center of gravity to point of contact with 
ground, measured perpendicular to reference line of 
fuselage, positive when contact is below reference 
line, ft (fig. l) 

angle of attack, relative to angle for zero moment, radians 
unless otherwise specified 


flight-path angle, radians 
generalized time, Cq_ 

ratio of ground drag load to ground vertical load (coeffi- 
cient of friction) 

ratio of M g to the angular inertia, Mgjiar^r 
mass density of air, slugs/cu ft 

angle between reference line of fuselage and ground, 
radians unless other-wise specified (fig. l) 
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Subscripts: 

a 

b 

e 

f 

max 

min 

o 


aerodynamic 

at first contact with the ground 
effective or equivalent 
fuselage reference line 
maximum 
minimum 

the point in time at which the impact of the tail is 
ended and the action of the rocker is begun 


CALCULATIONS 


Basis of the Calculations 

The calculations were concerned with the determination of loads and 
motions of an airplane landing and rocking on a curved undersurface such 
as that shown in figure 1. Most of the cases treated involved tail- low 
contact where the sequence of motion is of the type depicted in figure 2; 
here the spilling of wing lift provides an effective force for maintaining 
contact of the fuselage with the ground during damping of the pitching 
oscillations that ensue. In two cases the rocking was started by nose- 
low contact such as that shown in figure 3j here the initial increase of 
the wing lift requires careful shaping of the lower surface of the fuse- 
lage in order to maintain contact during stopping of the center of 
gravity, and, in fact, spilling of excess wing lift (natural or artificial) 
may be required to maintain ground contact after the contact point travels 
its maximum length. 

Equations used in the calculations are derived in appendix A. Equa- 
tions relating initial conditions for these equations to landing condi- 
tions prior to impulsive arrest of the point of contact are given in 
appendix B . In appendix C equations are developed which are approxima- 
tions to those in appendix A. These approximations were not used in the 
calculations but were useful in the analysis. The assumptions made in 
the calculations, for which equations are derived in the appendixes, are: 



1 


7 


impulsive arrest of the vertical motion of the effective mass at the 
point of initial contact, a rigid airframe, a rigid landing surface, 
constant coefficient of friction due to skidding, and no lag or effect 
of the ground on air lift or moments. Impulsive stopping of the con- 
tact point results in slight decrease of the center- of- gravity velocity 
before the contact point begins to travel along the length of the fuse- 
lage. A special shock absorber may be required to absorb the impulse 
in actual applications, but its design appears to be largely a matter 
of engineering and consideration of its details was not necessary in 
the analysis of subsequent loads and motions. 


Results of Calculations 

Tail-low landings .- Numerical results of the calculations for tail- 
low landings are summarized in table I. The data in this table are 
grouped according to four airplanes, for which the airplane character- 
istics are given in table II, shapes of the fuselage undersurface are 
given in figures 4 to and time histories of vertical acceleration, 
angular acceleration, and trim are given in figures 8 to 11. 

Most of the calculations were for a proposed bomber, which in some 
cases was considered to be equipped with an alternate wing which developed 
its full lift at higher angles. The second airplane, chosen for the pur- 
pose of obtaining an experimental check, was a fighter airplane for which 
a dynamic model was available. Calculations were made for original and 
modified lower- surf ace contours of this airplane. The third configuration 
was a transport airplane, the double-deck fuselage of which tended to 
increase the effect of ground friction and decrease vertical accelera- 
tions because of increased sweep- up of the rear of the fuselage. The 
fourth airplane was a research airplane for which the feasibility of 
fitting simple circular-arc rockers was studied. 

In those cases where table I refers to circular arcs or involutes 

(see fig. 12) these shapes were assumed and the accelerations and motions 

were calculated. The involutes (figs. 4(a) and 6) were very effective in 

reducing the high accelerations that occurred near zero trim for the 

circular-arc rockers (figs. 4(a), 7> 8(a), and ll)> tut gave high peaks 
at an intermediate point (figs. 8(a) and 10) and appeared to be desirable 
only in cases where the center region of the fuselage is weakened by bomb- 
bay doors. 

In those cases where table I refers to the type of load application 
the accelerations were assumed and the rocker shape calculated. Either 
the vertical center- of- gravity acceleration or the ground force was 
assumed to be constant, except that either the angular acceleration or 
the moment of the ground force about the center of gravity was not 
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allowed to exceed particular values in most of these calculations 
(figs. 4(b), 4(c) , 5, 8(b), 8(c), and 9). 

The landing conditions in table I represent cases in which there 
is no appreciable delay between stopping of the vertical motion of the 
effective mass at the tail and the beginning of forward motion of the 
point of ground contact. Provided the same conditions are maintained 
at the beginning of this forward motion, the histories of rocking can be 
interpreted in terms of landing conditions involving higher trims, 
higher values of wing lift, and either higher or lower sinking speeds, 
according to the initial value of the wing lift. Such interpretations 
involve the use of different starting points on histories such as that 
given in figure 13 for rotation about the tail point. Rotations of this 
type are calculated from the equations in the appendixes by substituting 
zero for the local radius of curvature of the rocker contour. Variation 
of the period of such rotation permits design of the contour for varia- 
tions of the landing trim. 

Nose-low landings .- Airplane characteristics and initial conditions 
used in the calculations for unflared nose-low landings are given in 
table III. Calculations were made for the actual contour of the air- 
plane and for the contour required to give a constant force of 15 percent 
of the weight. In the case of the actual contour, after the ground force 
dropped to zero at a time when the center of gravity had a sinking speed 
of 2 ft/sec, retraction of the flaps was assumed to occur in order to 
maintain contact. Similar spilling of the wing lift would be required in 
the case of the constant-force contour after the ground- contact point 
reached the tail. The two contours are shown in figure l4, and their 
histories of trim, center- of- gravity acceleration, and angular accelera- 
tion are given in figures 15 and 1 6 . 


EXPERIMENT 


Apparatus 

Tail- low landings were made with two configurations of a l/lO- scale 
dynamic model of the fighter airplane considered in the calculations. 

One configuration had the scale contour of the fighter, with the addi- 
tion of a l/k- by l/^-inch longitudinal strip that protected the under- 
surface of the model. The other configuration involved the addition of 
a variable- depth runner to obtain a modified rocker shape that had been 
calculated to give large reductions of both vertical and angular accel- 
erations. Some tests of this configuration were made with a l/b- by 
l/4-inch rubber strip insert to give some representation of the 
increased elasticity that might be expected for a full-scale airplane. 
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are plotted in figure 5 and shown in photographs of the model 
ures 17 and 18. 


in table ±, 
in fig- 


Landings were made by inclining a catapult to the proper flight. path 
angle and launching the model with the scale velocity. Time histories of 
normal acceleration, angular acceleration, and ground contact were trans- 
mitted to a fixed recorder by wires trailing behind the model. A dupli- 
cate run without these wires yielded no noticeable difference in the 
behavior of the model. The accelerometers had a natural frequency of 
100 cycles per second and were damped to 65 percent of critical. The 
ground- contact record was obtained by landing the model on metal plates 
placed in a concrete serviceway and using contact of a thin metal strip 
on the model to complete an electrical circuit. The metal plates were 
l/4-inch magnesium, except that after buckles 3/4 inch high were observed 
on a hot day the plates nearest the catapult were replaced with steel 
plates l/2 inch thick. 


Experimental Results 

Time-history records of center- of -gravity acceleration, angular 
acceleration, and ground contact are reproduced in figures 19, 20, and 
21 for landings of the previously discussed model configurations at a 
sinking speed of 9 ft/sec. The records for the modified contour were 
obtained with increased sensitivity of the instrumentation because lower 
accelerations were expected. Scale factors are indicated on the records 
and, in order to facilitate comparisons, the corresponding calculated 
histories and case numbers are also given. 

Notation is made on the records to indicate pulses which are due to 
dropping from the thick plate (full-scale drop of 2.5 inches) and due to 
joints between the plates. Such pulses cause degeneration of the rocking 
action into a bouncing action and, in the course of trial- and- error 
experimentation with the tail spring and its catch for preventing rebound, 
similar pulses frequently occurred following contact of the spring. The 
experiments, in general, confirmed the validity of the calculations, the 
expectation that careful selection of the tail shock absorber would be 
required, and the idea that irregularities of the landing surface should 
not be large in comparison with the combined deflections of the model and 
the landing surface. 

The path of the model indicated maximum lateral deviations of 
60 feet full scale, even though scaling of estimated cross winds to the 
full-scale airplane gave values up to 60 mph. In most cases there was 
a yawing of the model just before stopping, with maximum angles 


10 


approaching 90° for the basic configuration and about 45° for the modi- 
fied configuration. The stopping distance of the model indicated full- 
scale ground travel in the range of 1,000 to 1,200 feet. 


EFFECT OF MAJOR VARIABLES 


The analytical determination of the effect of major variables was 
complicated by the previously mentioned fact that the combination of 
reasonable shapes and desirable load histories was best achieved by a 
procedure in which the landing conditions were a result of the calcula- 
tion and could not be accurately controlled. Because of this complica- 
tion and the large number of variables involved, the following approxi- 
mate equations were derived in appendix B in order to focus attention 
on the major variables: 


• 2 / \2 
.. = ^b_Pe ] R 
J r \m / r 



These equations were derived by assuming that the conversion of sinking- 
speed energy into angular energy is the major action in rocking, and that 
the combined and partly compensating effects of air forces, weight, and 
skidding friction are negligible. Comments on some of these effects, 
with indication of conditions for which they become of major importance, 
are included in the discussions of the effect of y 2 /r, the significance 
of the scale radius of curvature R/r, and the relation between vertical 
and angular accelerations (determined by the variation of R/r with 
me/m) . In figure 22, the combination of these parameters in the approxi- 
mate equation for vertical acceleration is compared with previously pre- 
sented calculated acceleration histories for two sinking speeds. Repre- 
sentation of the major variations can be seen even though the differences 
due to the combined effects of air forces, weight, and ground friction 
are large enough to merit consideration. 

Effect of y 2 / r 

Since the parameter y 2 / r has the dimensions of acceleration one 
can expect that accelerations will be proportional to this parameter as 
shown by the approximate equations, and will be otherwise influenced by 
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various nondimensional ratios. These ratios Tna.v hp influenced, "by y— j~± 
ana thereby affect the proportionality indicated by the approximate 
equation, except that for a given vertical stroke in terms of r the 
average vertical acceleration must he proportional to y^/r. The 
buildup of angular velocity, the average angular acceleration, and the 
vertical and angular accelerations for a given point of contact are 
affected by the ratio of the total energy input to the sinking- speed 
energy and, according to the differences between these two energies, 
may deviate from proportionality to y2/r. 

Proportionality of maximum vertical acceleration to y^/r, for a 
given value of r, is shown in figure 23 for previously presented calcu- 
lations for sinking speeds of J>, 10, and 17 ft/sec. There is need for 
caution in assuming the general existence of such proportionality, as 
may be seen from the study in figure 24 of energy due to vertical motion 
of weight in excess of the wing lift. The scales in this figure deter- 
mine an equivalent sinking speed having the same amount of energy. Most 
of these equivalent sinking speeds are large compared with the design 
sinking speeds approximated with most of the data, thus indicating the 
importance of the subject energy relative to the sinking- speed energy. 

A semiempirical equation for determining this energy is shown in figure 24 
to be reasonably accurate . 

The great importance of energy due to deviations of the wing lift 
from the weight is further confirmed by the fact that when this energy 
subtracted from the sinking speed energy (nose- low impact), the net 
energy input into rocking, and the associated accelerations, were so low 
that difficulty was experienced in maintaining ground contact and the 
angular accelerations due to ground moments were less than those due to 
air moments (figs. 15 and 16) . In both nose-low and tail-low contacts if 
the sinking speeds are high enough the sinking- speed energy is predominant 
and the proportionality of acceleration to y2/r indicated by the approx- 
imate equation can be used without qualification. The sinking speeds 
required for this are less in the case where the energies add (tail- low 
contact), but on the basis of figure 24 they appear to be considerably 
higher than the design sinking speeds approximated in most of the 
calculations . 


Effect of R/r 

The approximate equations previously confirmed in figure 22 indicate 
that the accelerations are proportional to the ratio R/r. In this ratio 
the denominator r does not vary in a landing, but serves to establish 
the scale of the model. The numerator R locates an instantaneous hori- 
zontally translating center, directly above the point of contact, about 
which rotation occurs to give vertical acceleration at the point of con- 
tact which is simply the centrifugal acceleration t^R (see eqs. (3) 
and (43) ) • 
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For similar distributions of R and similar variations of r the 
maximum accelerations for a given r are proportional to R. In fig- 
ure 25 previously presented data for circular-arc contours show close 
agreement with this proportionality . Caution should he used in applying 
this result, however, since a change of R changes the trim range for a 
given travel of the contact point and changes energy sources which are 
affected ty trim. Previous study has indicated that large amounts of 
energy may be involved. 


Relation Between "t and y 

A minimum value of the maximum angular acceleration for a given 
buildup of angular energy in a given range of trim requires that ¥ be 
constant. The corresponding vertical acceleration approaches infinity as 
the moment arm of the ground reaction approaches zero. The condition for 
a minimum value of the maximum vertical acceleration with a given sinking 
speed and vertical travel of the center of gravity is a constant y. This, 
however, results in large angular accelerations for large moment arms about 
the center of gravity. For specified limits on both vertical and angular 
accelerations the sinking speed that can be tolerated will be a maximum 
if a period of constant ¥ for large moment arms is followed by a period 
of constant y for small moment arms . Such variations, and approximations 
thereto based on periods of constant ground moment and constant ground 
force, are included among the calculations previously presented. 

Maximum vertical and angular accelerations resulting from the various 
calculations for tail- low landings are compared in figure 2 6 through use 
of nondimensional acceleration coefficients derived in appendix C. In 
this plot of angular-acceleration coefficient against vertical-acceleration 
coefficient the solid slant line represents an approximate upper boundary 
associated with occurrence of maximum vertical and angular accelerations 
when the ground force has a maximum arm - for example, at the beginning of 
a case of constant vertical acceleration. At this time the aerodynamic 
forces and gravity are normally in approximate equilibrium and equa- 
tion ( 38 ) can be used to show that the slope of the limit line is equal 
to the maxim-urn value of (l + py)g/r or, with ground friction omitted 
from the approximate equations, Ig/r . This approximation indicates 
l/ r = 3*3 for "the limit line in figure 26. Actually the line was passed 
through the highest point, which is for a case of constant vertical accel- 
eration with l / r = 2.2. The 50-P ercerr t increase of slope from 2.2 to 3*3 
occurred primarily because the ground force due to deficiency of the ini- 
tial wing lift was 50 percent of the ground force associated with accel- 
eration of the center of gravity. 

The lower boundary of figure 26 approaches asymptotically a value 
determined by minimum angular acceleration associated with a constant 
value of this acceleration (expression (57))* Although this minimum is 
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reached only when the vertical acceleration is infinity, closeness of 
sume uf Ihe daLci Lu one ininiimim facilitated close approximation of the 
lower boundary. 

The left boundary in figure 2 6 is a curved dashed line fitted to the 
data and to the lower boundary. The combined boundary represents cases 
in which reductions in angular accelerations are achieved at the expense 
of increases in vertical acceleration, and vice versa. Although the mini- 
mum value of one of these accelerations can be achieved only at the 
expense of a large increase in the other, some of the data in the region 
of the dashed curve indicate that close- to-minimum values of both were 
achieved. Such cases appear to be optimum for design and, on the basis 
of data summarized in table I, offer acceptable values of both vertical 
and angular acceleration at design sinking speeds. 

Data for nose- low landings were omitted from figure 2 6 because there 
was no pronounced conversion of sinking- speed energy into angular energy, 
as was assumed in deriving the acceleration coefficients from the approxi- 
mate equations. These data (figs. 15 and 16) show a reduction of the 
angular velocity during stopping of the center of gravity after contact 
at the nose, and thus indicate the removal of energy from the airplane 
faster than the rate of reduction of sinking- speed energy. This result 
was caused by increases of wing lift which helped stop the center of 
gravity, by air moments, and by moments due to ground friction which 
opposed moments due to vertical ground force. The ability of these 
effects to absorb large quantities of energy is not surprising; large 
reductions of energy due to increase of wing lift might be expected in 
nose- low cases since in the opposite- rotating tail- low case the energy 
input due to decrease of wing lift beneath the weight was large in com- 
parison with the sinking- speed energy (fig. 24). Of less importance is 
the fact that energy input due to ground friction in a tail- low case 
caused a 20-percent increase of the maximum acceleration (fig. 10); 
reversal of this effect in the nose-low case accounts for an appreciable 
part of the difference between the two cases. 

The energy differences caused angular velocities and accelerations 
to be much less in nose-low landings than in tail- low landings. The 
higher sinking speeds that would be permissible would facilitate unflared 
nose-low contacts of airplanes which ordinarily would contact with their 
fuselage approximately level after a landing flare. The nose-low method 
is not feasible for airplanes which normally land at substantial trims, 
and it has the disadvantage that spilling of the wing lift may be 
required, as in the previously presented examples, to prevent the 
increased wing lift from causing loss of contact with the ground. The 
tail- low method is applicable to airplanes which normally land tail-low, 
and it has the advantage that decrease of the wing lift permits the 
weight to become effective in maintaining contact with the ground. 

Although angular velocities and accelerations are much greater in the 


i4 


tail- low case, they are acceptable even at sinking speeds typical of 
present design. These factors explain why the tail- low method has con- 
stituted the bulk of the present investigation. 


CONCLUSIONS 


1. On the basis of analysis and model tests, the method of landing 
and rocking on a fuselage undersurface of planned curvature to convert 
sinking-speed energy into angular energy for subsequent dissipation 
appears promising for special applications. 

2. In tail- low landings the accelerations are acceptable and the 
decrease of the wing lift caused by the rotation permits the weight to 
exert a powerful stabilizing force in maintaining ground contact and 
limiting angular rotation. 

3 . In tail- low landings the energy input due to descent of weight 
not supported by wing lift may be large in comparison with the sinking- 
speed energy. 

4. In nose-low contacts the increase of wing lift may be a predomir 
nant factor in stopping the vertical motion of the center of gravity, 
thereby greatly reducing angular velocities and accelerations, but 
spilling of the wing lift may be required in order to prevent rebound. 

5 . For a rigid fuselage landing on a rigid landing surface, the 
acceleration of the effective mass at the point of contact with the 
ground is simply the centrifugal acceleration t^R, where t is the 
angular velocity and R is the local radius of curvature of the 
undersurface . 

6 . In landings at sinking speeds high enough for ground forces to 
be predominant, minimum angular acceleration is achieved at the expense 
of large vertical acceleration, and vice versa, but it is possible to 
shape the undersurface so that close-to-minimum values of both are 
achieved. 

7 . The -contours required to give acceptable accelerations appear to 
be reasonable, but contours which otherwise appear to be reasonable will 
not give acceptable accelerations if there are large fluctuations of 
radius of curvature over lengths great enough for the fluctuation to 
integrate to differences in ordinate which are large compared with the 
local deformation of the fuselage and the landing surface . 

8 . The experimental data confirmed the expectation that irregulari- 
ties in the landing surface will be important if they are large in com- 
parison with deformations of the fuselage and landing surface. 
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9. The experimental data indicated that careful attention must be 
paid to the mechanism for stopping the vertical motion of the effective 
mass at the point of first contact so as to establish a smooth transi- 
tion to rocking on the curved lower surface. 

10. The experimental data showed satisfactory directional stability 
of a dynamic model during landing runout with initial rocking oscilla- 
tions induced by contact at a sinking speed which would induce limit 
loads in a landing gear. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., January 27, i 960 . 
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APPENDIX A 

EQUATIONS USED IN THE CALCULATIONS 


Rocker Dynamics 

The distance from the center of gravity to the ground plane (see 
fig. l) is given by the equation 


y = x sin r + z cos r 


( 1 ) 


Linear and angular velocities are related by the following equation, 
obtained by differentiating equation (l) and utilizing the geometric 
relation dz/dx = tan t: 


v = x cos r - z sin r (2) 

T 

where y/r is the distance from the center of gravity to the point of 
ground contact, measured along the ground. 

Linear and angular accelerations are related by the equation 

y = t 2 R - T 2 y + ( 3 ) 

T 


obtained by differentiation of equation (2) and use of the geometric 
relation 


R = & = - — ( 4 ) 

dT COS T T COS T 

Equation (3) corresponds to rotation of the airplane about an instan- 
taneous center located directly above the point of ground contact at a 
height equal to the radius of curvature R. The first term is the 
acceleration at the point of contact with the ground, which is simply 
the centrifugal acceleration about the instantaneous center. The second 
term corrects for the effect of vertical height in changing vertical 
components of the centrifugal acceleration about the instantaneous center. 
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The third term corrects the vertical acceleration of the center of 
gravity for the effect of angular acceleration about the instantaneous 
center . 

The force and acceleration normal to the plane of the ground are 
related by the equation 


F * my + mg - F a + F g 

The moment due to the vertical ground force on the rocker is 
M v = F(x cos t - z sin t) = F^ 

T 

and the moment due to ground friction is 


( 5 ) 


( 6 ) 


Md = kiyF 


(7) 


■where 

My moment due to vertical, ground force on the rocker 

Md moment due to drag load at the ground 

The angular acceleration is determined from the fact that the 
angular inertia moment is equal and opposite to the sum of the applied 
moments. Thus, 


-mr^T = M v + Md + Ma + Mg 


( 8 ) 


or, considering equations (6) and (j), 


F^yu + + Ma + M g 


mr 


( 9 ) 
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Aerodynamic Force and Moment 

The vertical aerodynamic force F a is given by the equation 

F a = L cos 7 + D sin 7 (lO) 


and the aerodynamic moment Mg^ is given by the equation 


M a = Moq + MqO, + M^.t 


(ID 


Nondimensional Equations 

The equations which have been presented are sufficient to permit 
solution by numerical integration. However, it is expedient to combine 
and simplify the equations and convert them to a nondimensional basis. 

Combining equations ( 5 ) , (9 ) > (lO), and (ll) results in the following 
equation: 


u 


iay + Fg + W - L cos 7 - D sin 7 


■)( 


yu + 


y\ 

il 

mr 2 


+ C. 


'“a' 


.aScdv 2 


+ + Cm. 


•+}Sc-£- £v 2 + m k 
2V 2 s 


( 12 ) 


Equation (12) can be simplified by assuming that (l) D sin 7 is 
negligible compared with L cos 7, ( 2 ) cos 7 = 1, ( 3 ) L = Cp^a&v 2 

where C]^ is a constant, (4) at first contact with the ground F a = W, 

( 5 ) moments are in balance at ground contact, (6) the pilot does not 
move the controls, ( 7 ) the flight-path angle 7 is represented by the 
equation 7 = y/V, and ( 8 ) the velocity V can be considered constant 
during a landing impact. 

With these assumptions, 

y - y b , . 

a = T - T b ^ — + % (13) 



d»0 


where the subscript b indicates values at the beginning of the landing 
impact . Equation (12) can then be reduced to the following form in 
which the terms are arranged so as to be nondimensional: 




• • \ 




7r _ l-y , F g , T " T b 

y - Yb\ U 

+ i_ ] + ^ 

T - T b 

y - h 

g \g w 

a b 

)\r 

vt) r 

a b 

%v 

+ (*-* + 

5Cm a' r 

5c ”i y 1 




A 2C “a 

a^yctbV 

2C % a b v2 

Wr 

i 




(15) 


A nondimensional equation of greater simplicity is obtained by 
introducing the generalized variables 


Y _ yy. 


T = rp 


(16) 

(17) 


and differentiating Y and T with respect to the generalized time 



(18) 


to obtain Y, Y, T, and T: 
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V 2 Ci 2 

Also, it is convenient to introduce the variables n and v where 

n - ^ 
lay 


2 - 

HITT 


With these substitutions, equation (15) reduces to 

(l - v)T = -Y(l + n) + T - Th - Cp^Y - Yb) ^Y + X 


where 


+ c 2 T - T b - C-l(y - Y b ) + C 5 T - C k Y 


Cl _ 

1 2 1 w/s / 


- dC™ 
Co = £. _S 
r dC L 


Cx = 


c m .) Cl 




In obtaining the equation of motion (15) and the simplified equa- 
tion (25), use was made of an equation for the applied force (eq. ( 5 ))* 
The first term in parentheses in equation (15) represents -F/W; that is, 
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* . # 


F = y _ £g 

w g " w 


T - T b + y - 


% 


a b V 


' (30) 


The first bracketed term in the simplified equation (25) represents 

-Cl 

-^10^; that is, 

= Y(1 + n) - T + Tb + C3 .(y - Y b ) 


(3D 


The motions and inertia loads resulting from the unknown rockers which 
give a particular history of ipa^ may be determined through use of 

equations (25) and (31) and a numerical integration procedure. Assumed 
values of the force must correspond to an upward push on the airplane 
since the ground is incapable of exerting a downward pull. 


The complete history of the motions and inertia loads of the 
unknown rockers which give a particular history of either Y or T can 
be obtained from equation (25) without use of equation ('31). However, 
when the accelerations are small or do not obviously correspond to an 
upward push of the ground on the airplane, equation (31) should be used 
as a check to see that the impossible condition of downward pull of the 
ground on the airplane is not implied by the assumed acceleration history. 


After determining the motion history of the unknown rockers which 
gives a desired load application, the unknown shapes may be obtained from 
the following equations, which come from equations (l) and (2): 


x _ Y 
r p 


T Y 

sin — + — cos 

P rp 


T 

M- 


(32) 


^Icosl-i 

r p p T 


T 

sin — 
M- 


(33) 


When loads and motions are to he determined for a rocker of known 
shape, the generalized vertical-acceleration parameter Y is determined 
by the generalized form of equation (3): 


y = p2 jj JL _y + p X 

pr T 


( 34 ) 
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This equation is used in combination with, equation (25) to obtain the 
complete loads and motions of a rocker of known contour. The radius R 
and the center- of- gravity height y may he explicit functions associated 
with a mathematical rocker or may he determined hy numerical means for a 
rocker of arbitrary shape . The ability of the rocker to maintain con- 
tinuous ground contact is determined by use of equation (31) to check on 
the sign of the computed force. In designs with small trim change during 
arrest of the tail, an approximate assumption of the force-time curve 
applied by the shock absorber will be sufficiently accurate for deter- 
mining r 0 . 

When a value of t q is determined a check should be made to insure 

that it is as large as the trim at which the rocker begins its working 
stroke. If this condition is not met the rocker will impact on a point 
closer to the center of gravity and receive a shock load. To allow for 
some variation in landing the rocker should be designed so that t q 
normally is greater than the trim at which the rocker cam becomes effec- 
tive . This introduces in the early part of the rocking a period of rota- 
tion about the end point of the rocker to which the equations for the 
cam- shaped rocker are applicable through the consideration of the radius 
of curvature at the point of ground contact as zero. 

The initial height of the center of gravity y Q may be determined 

from equation (l) by substituting the initial trim and the coordinates of 
the end point. This presupposes that the shock absorber is vertically 
positioned on the airplane so that the end point of the rocker contacts 
the ground at the completion of the arrest of this point. In order for 
a given shock absorber to fulfill this condition over a wide range of 
sinking speeds it should be of a constant- stroke type in which the 
applied forces are proportional to the square of the sinking speed. 

This force characteristic may be obtained from the vertically alined 
telescoping shock strut in which fluid is squeezed through an orifice 
at a rate proportional to the telescoping velocity. Such a shock strut 
also has the advantage that forces tending to extend the strut after 
arrest of the end point may be sufficiently small to prevent rebound and 
subsequent impact at another point on the rocker. This is particularly 
important when the trim after arrest of the tail is greater than the trim 
at which the rocker begins its working stroke, with the result that there 
is a period of rocking in which only aerodynamic moments due to motion 
set up by arrest of the tail are available for holding the tail on the 
ground . 
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APPENDIX B 
INITIAL CONDITIONS 


The initial conditions to be used in the equations of appendix A 
are those which exist after the vertical motion of the tail has been 
stopped by the tail shock absorber. These conditions (indicated by 
the subscript o) are related to those at the beginning of tail impact 
(indicated by the subscript b) by momentum considerations related to 
the impulse required to stop the vertical velocity of the tail at 
contact with the ground. This impulse is equal to ir^y^, where mg 
is the effective mass at the tail. The change in vertical momentum 
of the center of gravity is also equal to this impulse, so that 

m(y b - y Q ) = m e y t ^ (35) 

or 

y Q = y b - ^t,b ( 56 ) 

where y b ^ b may be obtained from the equation 

y t = y - Ti 

Here l is the horizontal distance from contact point to center of 
gravity and me/m is determined as follows. 

The inertia of the effective mass at the tail and the center of 
gravity are both equal to the ground force and are equal to each other. 
Thus, 


®e*yt = m y (57) 

where 

y t = y - Ti 

The angular acceleration is related to the vertical acceleration hy 
the condition that the sum of the angular- inertia moment, the vertical- 
force moment, and the ground friction moment is equal to zero. Thus 


2k 


or 


Then 


From equation (37 ) > 


r^mr 2 ^ + F(Z + |ay) = 0 
T^mr 2 ) + my( l + |ay) = 0 


•f = _ lil i t*y) 

■v.2 



(38) 


Substituting equation ( 38 ) gives 

me = 1 

m “ 1 + (39) 

or, when q. during the tail impulse is equal to p during the subse- 
quent rocking. 


me 

m 


m 


1 + 


vT, 


+ 4)y 


*/ 


The trim r n is given by the equation 


T o - T b + Vi + 


t dt dt 


(41) 


where the duration of the tail impulse tq and the integral are deter- 
mined by the design of the tail shock absorber. 
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APPENDIX C 

APPROXIMATE VARIATION OF VERTICAL 
AND ANGULAR ACCELERATIONS 


The variation of the vertical and angular accelerations can he 
approximated by assuming that the gain in angular energy is equal to 
the loss of sinking- speed energy (negligible or compensating effects 
of unbalanced weight, skidding friction, and air moment) . Such an 
assumption is most reasonable for design impacts which involve maximum 
sinking-speed energy and maximum accelerations, and hence are of the 
greatest interest. With the above assumption, 


2*2 

mr^r 



and, with y = tZ, 


2*2 * 2 

= m e y^ 




where (see eq. (59)) 


“e = 



The vertical acceleration at any instant may be determined on the 
basis that the instantaneous center of rotation is directly above the 
point of ground contact at a distance equal to the radius of curvature 
of the rocker at this point (eq. (5))* Thus the vertical acceleration 
of the effective mass is simply the centrifugal acceleration due to 
this radius - the term t^R of equation (5) • The acceleration at the 
center of gravity is related to the acceleration of the effective mass 
through the fact that both inertia forces are equal since they are both 
equal to the vertical force. Therefore 


•• • 2-d 

my = n^T R 


(^ 3 ) 
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Combining equations (42) and (43) gives 



(44) 


The angular acceleration is determined by the relation that the 
angular inertia is equal to the applied torque. This may be expressed 
as follows, where the inertia of the effective mass has been substituted 
for the ground force: 


mr^r = mg-f^Rl 


( 45 ) 


Combining equations (43), (44), and (45) yields the following equation 
for angular acceleration: 



vl/2 

R 

m / 


(**) 


These approximate equations permit the calculation of acceleration 
histories without the necessity of performing the numerical integrations 
required "by the more exact equations used in the present investigation. 
They show that for a given point of contact the radius of curvature of 
the lower surface contour at that point is the sole geometric parameter 
affecting the accelerations, and that accelerations are proportional to 
this radius of curvature. Decrease of the radius of curvature may he used 
to reduce the associated ground forces to any desired value, hut excessive 
use of this means will result in angles of the contour slope greater than 
the landing trim and thereby cause contact to occur at a point which is 
closer to the center of gravity. This increases the effective mass at the 
point of initial contact, and the increased shock of initial contact tends 
to offset the subsequent reduction in force. To reduce the shock of 
initial contact to a minimum and insure contact at the same point over a 
reasonable range of landing trims so that a special shock absorber can be 
effective, it will be generally desirable to limit the integrated effect 
of curvature so that 



^ T min 
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I’d' a given strength distribution of the airplanp. the optimum 
lower surface contour will generally he that for which the curvature 
integrates to the indicated maximum and is so distributed that the 
critical sinking speed causes a critical load to occur at some point at 
al 1 times. In comparing results for different curvature distributions, 
corrections should be made for the effects of integration to different 
maximum slopes. With a view toward this purpose let 



(47) 


where 


H c 



characteristic radius of curvature 


function describing the curvature distribution 


Equation ( 4 ) can be reduced to 


dx 

dr 


(48) 


on the "basis that the range of cos t is very close to unity. Combining 
equations (Vf) and (48) yields 


dr 



(49) 


from which it is apparent that 


r 



(50) 


Combining this with equations (44) and (46) and considering that, for a 
particular curvature distribution, R at a particular x/r is propor- 
tional to R c , gives 


yT 


oc 



r 


(51) 
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y-u 

TrT « (52) 

r 

The proportionality of acceleration to y^ 2 indicated in 

equation (52) stems from equation (42) which assumes that the angular 
energy is equal to the loss of sinking- speed energy. Since a semiem- 
pirical expression for energy due to unbalanced weight has been pre- 
sented (fig. 24), use of it to modify expressions (51) and (52) appears 
desirable. From equation (42), 

| mr2 'to = |myb 2 (55) 

and it is desirable now to use the equation 

^ f max = l^b 2 + E w (5*0 

The desired change will be effected if y^ 2 in equation (55) is 
replaced by y^ 2 + 2g-~. Using the equation for E w from figure 24 and 

rearranging expressions (51) and (52) to isolate constants of propor- 
tionality results in the following acceleration coefficients representing 
these constants: 


y^ 2 + 2gl 


J max 

0.7^ - 

w 


rr 


0.5 


b'min [ 
W / 



(55) 


w r2T 


y* + 2gll 


0 . 7 ^- 0.5 


^ r min\ 


W 


■j(yo ~ ymin 


(56) 


These acceleration coefficients are plotted against each other in fig- 
ure 26, which presents the results of various calculations that were made. 

Through use of equation (54) the angular-acceleration coeffi- 
cient (5 6) can be written 
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T max T 


- O 

T*~ 


max 


(57) 


From this form of the angular-acceleration coefficient one can determine 
that its minimum value occurs for constant angular acceleration and is 
equal to l/2. Thus 


• • Q 

max 1 T 


y b 2 + 2 g(l - 0.7^ - 0.5^)(y o - y^ n ) 


>1 
“ 2 


(58) 


(within the limits of the various assumptions that have been made) 
defines the lower boundary in figure 26 if r is converted to degrees as 
in this figure. 
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TABLE I 


va- 


^ rtTTT A rrrr Awn TiAT* 


LANUIIJC C 


Case 

Sinking 

speed, 

h> 

ft/sec 

Initial 
trim, a 
T b> 
deg 



Initial 

hi* 

Coefficient 
of friction, 
4 

Rocker type 

Maximum 

vertical 

acceleration, 

y/g 

Maximum 

angular 

acceleration, 

radians/sec^ 

Maximum 

negative 

trim, 

T, 

deg 

Bomber airplane 

1 

3 

7 

1 

0.4 

Involute 

0.25 

- 0 . 25 , 0.43 

3.5 

2 

10 

7 

1 

.4 

Involute 

1.37 

-.98, 1.27 

5.8 

5 

17 

7 

1 

.4 

Involute 

4 .oo 

- 2 . 4 , J.J 

8.0 

4 

17 

8.6 

1 

.4 

Circular arc 

5.6 

" 2.3 

— 

5 

9.2 

5.6 

1 

.4 

Limited moment* 5 

•9 

-2.0 



' 6 

19-8 

16.5 

•9 

.4 

Limited moment* 5 

1.5 

- 1.65 

— 

7 

17-7 

12.5 

.8 

.4 

Limited moment* 5 

3.0 

-1.8 

— 

8 

25-8 

19.6 

1 

.4 

Limited moment* 5 

6.0 

-2.83 

— 

9 

10.1 

6.2 

1 

.4 

Constant acceleration 

•9 

" 3-6 



10 

13.2 

9-7 

.6 

.4 

Constant acceleration 

•9 

- 3.2 



11 

13.1 

l 4 . 2 

• 7 

.4 

Constant acceleration 

.67 

- 3.3 



12 

29.O 

7 

• 7 

.4 

Constant acceleration 

6.0 

-17.1 

— 

Fighter airplane 

13 

9 

13 

l 

0.4 

Unmodified 0 

io .4 

-16, 21 

10.7 

14 

9.6 

12.5 

l 

.4 

Limited moment 0 ’ d > e 

2.5 

-6.7, 12 


15 

9.6 

12.5 

l 

.4 

Limited moment 0 

2.5 

-6.7, 6.7 

— 

16 

9.6 

12.5 

l 

.4 

Limited moment'*’ ** 

2.5 

-6.7, 24.6 

— 

17 

9.4 

11.4 

.9 

.4 

Limited moment^ 

3-3 

" 6.7 

— 

Transport airplane 

18 

10 

13 

1 

0.4 

Involute, R c = 0 ft 

1.0 

-1.07 



19 

10 

15*5 

1 

.4 

Involute, R c = 33 ft 

1.13 

-1.0, 1.52 

9 

20 

10 

14.5 

1 

.4 

Involute, R c = 120 ft 

1.35 

-1.0 

— 

21 

10 



14.5 

1 

0 

Involute, R c = 120 ft 

1.0 

-.77 

— 

Research airplane 

22 

10 

11.8 

1 

0 

Circular arc, 88- ft radius 

8.1 

- 15.6 


23 

10 

17.2 

1 

0 

Circular arc, 57 " ft radius 

5.6 

-10.9, 11.4 

13 


a Based on 1 ° change in trim between tail contact and start of forward movement of contact point. 
^Moment of the vertical ground force limited to that produced by W /2 applied at the tail. 

0 Checked by experiment. 

= y*f when t is less than indicated value. 
e Flaps assumed to be raised at zero trim. 

f 

Forebody computed without limitation of moment. 

^F = 4 W when t is less than indicated value. 






TABLE II 


AIRPLANE CHARACTERISTICS USED IN COMPUTATIONS 
FOR TAIL- LOW LANDINGS 




Bomber 

Fighter 

Transport 

Research 

W/S, lb/ sq ft . . . 

« f • 

7^ 

46 

78 

(a) 

V, knots 

• • • 

112 

100 

112 

(a) 

°L a 

• • • 

• 5 . 1 , b 3.2 

4.2 

4.1 

° 3 C L 

Of for L = W, deg . 


9 , b 25 

l 6 

16 

14 , C 19 

h> ft 

• • • 

55 

18.6 

70 

16.9 

r, ft 

• • • 

25 

6.6 

22 

5-9 

c, ft 

• * • 

11 

8.09 

21 

(a) 

C m a , 

• ♦ • 

- 0.4 

-0.3 

-0.4 

(a) 

Cm*! • * 

• • • 

d -40 

-7.14 

d -l 8 

(a) 

Cm 

. . . 

(e) 

-2.86 

(e) 

(a) 


^his characteristic was removed from consideration by assuming 
C T to have the indicated relation to the Ct at landing, by assuming 

changes in 7 to be small compared with changes in t, and by omitting 
changes in aerodynamic moment. 

^ Alternate wing. 

Alternate wing or lower landing speed, 
d 


'Assumed to be of negligible importance. 
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TABLE III 

AIRPLANE CHARACTERISTICS USED IN CALCULATIONS 


FOR NOSE- LOW LANDING 


jRaising of flaps assumed to have negligible effect on trimming 
moments, or to have its effect canceled by control movement.! 


W, lb . . 
V, knots 
y, ft /sec 
S, sq ft 
c, ft . . 
r, ft . . 


^Cm 

dC L 

C„ 

“a 

C m| 


58,750 

lkO 

9 

1,175 

16.8 

10.9 

4.5 


-0.13 


-0.58 

-2.9 

- 7.1 



Figure 2.- Motion sequence in a tail- low rocker landing 
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Figure 6.- Rockers for transport airplane. 




Trim, % f deg Angular acceleration, radiana/sec Vertical acceleration at 





e 

al 
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Time, t, sec 


(c) Constant acceleration. 
Figure 8.- Concluded. 






Time, t, sec 


Figure 10.- Time histories of vertical acceleration, angular accelera- 
tion, and trim for transport airplane with a sinking speed of 10 feet 
per second. 





Figure 13 
tail pi 
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Figure 15*- Load and motion histories for nose-low landing. 

Unmodified contour. 



7» g units 


y, « 

7> g units 
T, radians/sea 
T, radians/sec 2 
T, deg 
7, ft/snc 



-8 I t I I 1 1 J 1 

0 .08 .16 .24 . 32 . 40 .48 .56 

Time, t, sec 


Figure 16.- Load and motion histories for nose-low landing 

Modified contour. 
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L- 86292 

Figure 17 l/lO- scale dynamic model of a fighter with a rutting strip fitted 

to the fuselage contour (case 13)* 



L-90408 

Figure 18.- l/lO-scale dynamic model of a fighter with a computed rocker (case l4) . 




Ground contact 


✓ • 



Figure 19 •- Ground contact, vertical acceleration, and angular acceleration for landing 
of a model of a fighter on a rubbing strip fitted to the lower surface of the fuselage. 
Case 13 in table I. 












Approximate equation: 








Figure 23.- Variation of vertical acceleration with sinking speed. 



Figure 24.- Evaluation of approximate equation for the effect of unbalanced weight 








